Mott transition in the Hubbard model on the anisotropic kagome lattice 
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We investigate the Mott transition in the anisotropic kagome lattice Hubbard model using the 
cellular dynamical mean field theory combined with continuous-time quantum Monte Carlo simula- 
tions. By calculating the double occupancy and the density of states, we determine the interaction 
strength of the first-order Mott transition and show that it becomes small as the anisotropy increases. 
We also calculate the spin correlation functions and the single-particle spectrum, and reveal that 
the quasiparticle and magnetic properties change dramatically around the Mott transition; the spin 
correlations are strongly enhanced and the quasiparticle bands are deformed. We conclude that 
such dramatic changes are due to the enhancement of anisotropy associated with the relaxation of 
frustration around the Mott transition. 

PACS numbers: 71.30.+h 71.10.Fd 71.27.+a 



Geometrically frustrated electron systems have pro- 
vided hot topics in the field of strongly correlated elec- 
tron systems, and have uncovered various new aspects of 
the Mott transition. The discovery of heavy fermion be- 
havior in LiV204 P, Q with pyrochlore lattice structure 
has activated theoretical studies of electron correlations 
with geometrical frustration [3-6]. More remarkably, re- 
cent experiments on the triangular lattice organic mate- 
rials /^-(BEDT-TTF)2X have revealed a novel spin liquid 
ground state in the Mott insulating phase |7l-[ll|. 

The kagome lattice (see Fig. [1] (a)) is another pro- 
totype of frustrated systems, which shares some essen- 
tial properties with other frustrated lattices. The lo- 
calized electron systems on this lattice have been inten- 
sively studied and many unusual properties have been 
found [l2|. In particular, a spin liquid state observed 
in the herbertsmithite ZnCu3(OH)6Cl2 with kagome lat- 
tice structure [l3|, Q has received considerable atten- 
tion. A related spinel-oxide Na4lr308 with hyperkagome 
structure, which is a three-dimensional analog of kagome 
lattice, is also proposed as a candidate for the spin liq- 
uid (il, [iH. The issue of electron correlations on the 
kagome lattice was addressed by applying the FLEX ap- 
proximation [17] and the QMC method [18] to the Hub- 
bard model. Also in our recent study, we found a heavy 
Fermi liquid state emerging in a metallic phase near the 
Mott transition [19]. Furthermore, an idea of chirality- 
spin separation was proposed to explain low-energy char- 
acteristic properties [20|, HH . In spite of such intensive 
theoretical investigations, however, the effects of spatial 
anisotropy have not yet been addressed. The system- 
atic exploration of the anisotropy is desired to elucidate 
how the frustration affects the nature of Mott transition, 
the quasiparticle formation, the magnetic properties, etc. 
It has indeed been clarified that in a moderately frus- 
trated system on the anisotropic triangular lattice, the 
Mott transition shows a reentrant behavior [qI, HH, |23|, 
where the low-temperature transition is governed by the 
enhanced antiferromagnetic (AFM) correlations. This is 
a unique feature that may occur neither in the unfrus- 



FIG. 1: (a) Anisotropic kagome lattice, (b) effective cluster 
model used for the CDMFT and (c) first Brillouin zone of the 
anisotropic kagome lattice (c). 



trated square lattice [2J nor in the fully frustrated tri- 
angular lattice [ill . This naturally motivates us to study 
the effects of frustration on the Mott transition by sys- 
tematically controlling the spatial anisotropy. 

In this paper, we study the Mott transition in the Hub- 
bard model on the anisotropic kagome lattice by means 
of the cellular dynamical mean field theory (CDMFT) 
[26| combined with the continuous-time quantum Monte 
Carlo (CT-QMC) method [27]. We first treat the 
isotropic system to confirm our previous results studied 
with a slightly different method [19]. As the anisotropy 
is introduced, electronic properties drastically change 
around the Mott transition; the spin correlations are 
strongly enhanced and the quasiparticle bands are de- 
formed substantially even if the system is slightly away 
from the isotropic point. We elucidate that such dra- 
matic changes are due to the enhancement of anisotropy 
driven by the relaxation of frustration around the Mott 
transition. 

We consider the Hubbard model on the anisotropic 
kagome lattice shown in Figdja). Here, we set the lattice 
geometry as a square form, where the first Brillouin zone 
is — TT < kx^ky < TT diS shown in FigHJc). The Hamilto- 
nian reads 
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FIG. 2: Double occupancy as a function of interaction 
strength U for several at T — 0.05. 

where {cja) creates (annihilates) an electron with spin 
cr at site i and n^cr = c\(j^icr- The hopping integrals and 
the Hubbard interaction are denoted by t {t') and [/, 
respectively. The system corresponds to the fully frus- 
trated kagome lattice at t' /t = 1, and frustration be- 
comes weaker with decreasing t' /t. The end member at 
t' /t = is called a decorated square lattice. We choose 
the energy unit as t = 1 hereafter. To investigate strong 
correlations and geometrical effects, we use CDMFT [26], 
a cluster extension of DMFT [28]. CDMFT has been suc- 
cessfully applied to a lot of frustrated electron systems 
so far [ii,|22|,[2i,|2i,|30j. 

In CDMFT, the original lattice problem is mapped 
onto an effective cluster problem. Each unit cell of the 
anisotropic kagome lattice has three sites, as shown in 
Fig. m^a). We thus end up with a three-site cluster 
model coupled to a self-consistently determined medium 
illustrated in Fig. [IJb). Given the Green's function for 
the effective medium, ^, we compute the cluster Green's 
function G and the cluster self-energy S in the effective 
cluster model, where O denotes a 3 x 3 matrix. To cal- 
culate G and S in the effective cluster model, we use the 
CT-QMC method. In CT-QMC simulations for an ef- 
fective cluster model, the sign problem strongly depends 
on the choice of basis set in the local Hamiltonian. We 
choose the basis set which diagonalizes the local hopping 
matrix tioc written in the sublattice basis, 

( {) t t\ 

hoc = \ t ^ t' \. (2) 

\t t' {) I 

We find that this choice substantially reduces the sign 
problem in QMC simulations. We iterate the DMFT self- 
consistent loop until the convergence of this procedure is 
achieved within 50 iterations at most. In each iteration, 
we typically use 2.5 x 10^ QMC sweeps to reach sufficient 
computational accuracy at very low temperature, T = 
0.05. 

Let us now investigate the Mott transition on the 
anisotropic kagome lattice at half filling. In order 
to find evidence of the Mott transition, we calculate 
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FIG. 3: Density of states for typical U and at T = 0.05. 



the U dependence of the double occupancy Docc. = 
I Y^m=i {"^rn^^mi) for various t' . The double occupancy 
monotonically decreases with increasing U and shows a 
discontinuous jump at the interaction Uc We also find 
hysteresis, which signals the emergence of the first or- 
der Mott transition. In the isotropic kagome system for 
t' = 1.0, the Mott transition occurs at fairly large in- 
teraction strength Uc ~ 8.4 compared with the band- 
width ly = 6. This is consistent with our previous 
study using CDMFT combined with Hirsch-Fye QMC 
[19]. In anisotropic cases, the interaction strength Uc be- 
comes small as t' decreases. Note that the band width 
W/t = 4^/2 ~ 5.66 for = and W is almost unchanged 
when t' varies. Therefore, we can say that the metallic 
region shrinks as geometrical frustration is weakened. 

To see the metal-insulator transition more clearly, 
we calculate the density of states (DOS) A{uj) = 

-^YL.=i^'^^rnm{^ + i^) by applying the maximum 
entropy method (MEM) [31] to the imaginary-time QMC 
data. In Fig El we show the DOS for t' = 1.0, 0.8 and 0.6. 
We clearly see that the insulating gap opens around the 
Fermi level in the large U region beyond Uc For t' = 1.0, 
the three bands are strongly renormalized and the heavy 
quasiparticle states emerge near the Fermi level in the 
strong correlation regime. A similar peak structure also 
appears for = 0.8, 0.6 in the weaker U region than that 
for t' = 1.0. 

Let us now turn to the magnetic correlations. To see 
how frustration affects them around the Mott transi- 
tion, we calculate the nearest-neighbor spin correlation 
function (5'f ^f+i). Here, i denotes a sublattice index as 
shown in Fig. [TJ In Fig. 21 we show {S^S^_^i) as a func- 
tion of U for several t' . In the isotropic case, {Si 82)^ 
(SfS^) and (/S'l^S'l) are equivalent and have small nega- 
tive values for any U. This indicates that the AFM spin 
correlation is suppressed due to the strong frustration. 
As finite anisotropy is introduced {f = 0.9,0.8,0.7,0.6), 
the spin correlations dramatically change near the Mott 
transition; {SIS2) and {SIS^) show strong AF correla- 
tions while (S'I'S'I) tends to be ferromagnetic. These be- 
haviors indicate the tendency to a ferrimagnetic ordering 
(or an AFM- type ordering on a decorated square lattice). 
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a function of U in the metallic phase. 



FIG. 4: Nearest-neighbor spin correlation function {SfSfj^i) 
as a function of U for several t\ where the suffix i specifies a 
sublattice as shown in Fig. [T] 
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FIG. 5: Momentum resolved spectral weight at the Fermi 
level. In (c) and (f), we show the Fermi surfaces deduced 
from (b) and (e) together with those for the noninteracting 
case for reference. 



We thus find that magnetic properties in the system are 
very sensitive to the anisotropy, whose effect is partic- 
ularly enhanced in the vicinity of the Mott transition 
point. 

The relaxation of frustration strongly affects the for- 
mation of quasiparticle bands. We next investigate 
the momentum resolved spectral weight at the Fermi 
level, A{\<i,UJ = 0) = -^YL.=i^'^^rnm{^,^ 



0), in the metallic phase. Here, we approxi- 
mately compute A(k, a; = 0) as A{\<i^uj = 0) ~ 

-^YL.=l^'^^rnm{^,'l^n ^0). In Fig. O WC shoW 

color plots of 74(k, = 0). Also, shown in Fig. [5] (c) 
and (f) is the Fermi surface deduced from the trajectory 
of their peaks. In the isotropic case, A(k, = 0) for 
U = 4.2 has strong intensity at the Fermi surface and 
the profile of the Fermi surface is consistent with that 
in the noninteracting case. Near the Mott transition, for 
U = 8.3, the Fermi surface keeps the same shape as in 
the noninteracting case, although the quasiparticle peak 
in 74(k, cj) is somewhat obscured at finite temperatures. 
In contrast, for t' = 0.8, the geometry of the Fermi sur- 
face drastically changes near the Mott transition. In the 
weak coupling regime, for U = 3.65, the Fermi surface 
has a shape similar to that in the noninteracting case. As 
U increases, for U = 7.2, the shape becomes elongated 
and quite different from that in the noninteracting case. 
It is known that the shape of the noninteracting Fermi 
surface is elongated along the kx = ky direction as t' de- 
creases [17]. Therefore, we find that the band structure 
of quasiparticles becomes much more anisotropic due to 
the enhancement of hopping anisotropy, which is induced 
by the relaxation of frustration. 

In order to confirm the above scenario, we investigate 
the renormalization of hopping integrals. In the low- 
energy limit, the Green's function in the effective cluster 
model may be given by 

G{lo) = Z'/' [{lo + - tl, - Z^/^ (3) 

where t^^^ is the matrix representation of the renor- 
malized hopping integrals defined by t*^^ = Z^/^(tioc — 
Ret{uj = 0))Zi/2, where Z = {1 - dRet{uj)/duj\^=o)~^ • 
In tJ^Q^, the original hopping integral t {t') is renormalized 
as, t ^ {t' t'^). The renormalized hybridization 
function is defined similarly: A*(co') = Z^I^K(uS)Z^I^ . 

In Fig. [6l we show the /7-dependence of the ratio 
of the renormalized hopping integrals t"" /f". In the 
isotropic case, is always unity due to the symme- 

try requirement, giving rise to strong frustration near the 
Mott transition. For the anisotropic cases of t' /t = 0.9 
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and 0.8, does not show any drastic change in the 

small U region. However, in the vicinity of the Mott 
transition, rapidly decreases, implying that the 

anisotropy in hopping integrals is considerably enhanced 
there. As a result, geometrical frustration becomes weak. 
In other words, strong frustration developed in the nearly 
isotropic models triggers the strong renormalization of 
anisotropic hopping in the vicinity of the Mott transi- 
tion. The enhanced anisotropy deforms the band struc- 
tures and develops the spin correlations. We thus con- 
clude that the dramatic changes found around the Mott 
transition are due to the enhancement of anisotropy as- 
sociated with the relaxation of frustration. 

We note that the above mechanism for relaxation of 
frustration may be generic for frustrated electron sys- 
tems in the vicinity of the Mott transition. For example, 
the reentrant Mott transition found for the anisotropic 
triangular lattice compound [9] may be caused by the 
above mechanism. Namely, if the lattice anisotropy is 
weak, a metallic state near the Mott transition keeps 
the isotropic nature approximately in the weak coupling 
regime (high temperature) and thus behaves as a fully- 
frustrated model. However, as the temperature is low- 
ered, the system enters the strong-coupling regime, where 
the anisotropy is enhanced and the AFM correlations are 



developed. This may trigger the reentrant Mott transi- 
tion where the low-temperature Mott phase is accompa- 
nied by enhanced AFM correlations. 

In summary, we have studied the Mott transition 
in the Hubbard model on the anisotropic kagome lat- 
tice by means of the CDMFT combined with CT-QMC. 
For anisotropic lattice systems, we have found that the 
quasiparticle and magnetic properties drastically change 
around the Mott transition; the spin correlations get 
strongly enhanced and the quasiparticle bands are de- 
formed. It has been elucidated that these behaviors are 
due to the relaxation of frustration caused by the en- 
hancement of anisotropy around the Mott transition. 
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